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摘要 ”为 更 好 的 理解 和 掌握 熔 盐 换 热 器 的 运行 特征 ， 并 给 炊 盐 堆 中 换 热 器 的 设计 和 运行 累积 经 验 ， 本 文 以 让 
国 橡树 岭 国家 实验 室 10MWt 熔 盐 实验 堆 (MSRE) 用 的 主 换 热 器 为 研究 对 象 ， 基 于 该 换 热 器 的 设计 参数 ， 应 
用 管 壳 式 换 热 器 的 理论 计算 方法 (包括 Kem 方法 和 Bell-Delaware 方法 )、 换 热 器 设计 软件 HTRI Xchanger Suite 
和 计算 流体 力学 方法 (Computational Fluid Dynamics, CFD) 对 换 热 器 的 关键 性 能 指标 〈 如 换 热 系数 、 压 降 和 
换 热 功率 ) 分 别 进行 了 计算 分 析 ， 并 与 MSRE 主 换 热 器 的 实验 数据 进行 比较 分 析 。 结 果 表 明 : 两 种 理论 计算 
方法 (Kem 方法 、Bell-Delaware 方法 )、HTRI 软件 及 CFD 仿真 模拟 的 结果 与 实验 结果 的 差异 均 在 可 接受 范 目 
内 ， 其 中 Kem 方法 与 MSRE 实验 所 得 的 换 热量 相差 最 大 、 约 为 15%，HTRI 软件 模拟 与 MSRE 实验 所 得 的 总 
换 热 系数 相差 最 小 、 为 0.16%; 理论 计算 方法 、 换 热 器 设计 软件 及 CFD 仿真 均 不 同 程度 适用 于 熔 盐 - 熔 盐 换 热 
器 设计 ， 且 CFD 仿真 还 可 对 换 热 器 内 部 传 热 和 流动 细节 进行 直观 描述 。 
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Simulation analysis of heat transfer and flow characteristics of a U-tube heat exchanger in a 


molten salt reactor 
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Abstract [Background]: The primary heat exchanger (PHX) used in the 10M Wt Molten Salt 
Reactor Experiment (MSRE) at Oak Ridge National Laboratory (ORNL), is a U-tube heat exchanger, 
where the shell side (hot side) contains the fuel salt from the primary loop, while the tube side (cold 
side) carries the coolant salt from the secondary loop. [Purpose]: This study aims to deepen the 
understanding and mastery of the operational characteristics of molten salt heat exchangers, and to 
accumulate experience in their design and operation within molten salt reactors. [Methods]: the 
MSRE-PHX is modeled based on the design parameters, theoretical calculations for shell and tube 
hear exchanger (Kern method and Bell-Delaware method), software simulation (HTRI Xchanger 
Suite) and computational fluid dynamics (CFD) simulation are performed, critical performance 
metrics, such as the heat transfer coefficient, the pressure drop, and the heat transfer power, are 
obtained and compared to the MSRE operation data. [Results]: The findings indicate that the 
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discrepancies from theoretical calculations, HTRI software, and CFD simulations, are all within 


acceptable margins to the experimental data. Notably, the greatest variance was found with the Kern 


method, which showed a deviation in heat transfer quantity of about 15%, while the smallest 


discrepancy was observed in the overall heat transfer coefficient calculated using HTRI software, 


differing by merely 0.16% from the experimental data; [Conclusions]: All of the methods are 


suitable and applicable for designing and studying a molten salt shell and tube heat exchanger; 


moreover, the CFD simulation can provide fine localized details of the heat transfer and flow of the 


molten salt fluid. This offers substantial theoretical support and practical guidance for the future 


design and improvement of molten salt heat exchangers. 
Key words Molten Salt Heat Exchanger, Kern Method, Bell-Delaware Method, HTRI Xchanger 


Suite, CFD Simulation 


随 着 全 球 能 源 需 求 的 不 断 增 长 和 环境 问题 的 
国际 核能 界 的 共识 。 在 第 四 代 反 应 # 


盐 堆 可 在 700°C 以 上 的 高 温 条 作 
远 高 于 传统 压 水 堆 33% 左 右 的 效率 水 习 


T Tis e FR 


液态 燃料 熔 盐 堆 中 ， 含 有 裂变 材料 的 燃料 盐 在 一 回路 9 
PHX) 将 热量 传递 给 二 回路 冷却 盐 。 
流体 力学 问题 。 因 此 ， 对 炊 盐 换 热 器 中 流动 与 传 热 特性 的 研 


出 ， 开 发 安全 、 高 效 、 清 洁 的 新 型 核能 系统 已 成 为 


益 突 
全 中 ， 熔 盐 堆 因 其 固有 的 安全 


性 、 高 热效率 和 多 功能 性 备 受 关注 让。 熔 


这 一 过 程 不 仅 涉及 高 温 


系统 中 的 关键 设备 之 一 ， 其 性 能 直接 影响 着 熔 盐 堆 的 安全 性 、 经 济 性 和 可 靠 性 所 。 在 
循环 流动 ， 通 过 主 换 热 器 (Primary Heat Exchanger, 
下 的 材料 和 结构 安全 ， 还 涉及 复杂 的 热力 学 和 
f 究 至 关 重 要 。 


F 下 运行 ,实现 高 效 紧 凑 的 布雷 顿 循环 发 电 ， 其 理论 发 电 效 率 有 望 达到 44%， 


熔 盐 的 热 物 性 ， 包 括 密 度 、 粘 度 、 比 热 容 和 热 导 率 等 ， 是 换 热 器 设计 的 关键 参数 。 在 20 世纪 的 熔 盐 堆 


实验 (Molten Salt Reactor Experiment, MSRE) 等 项 目 中 ， 美 


国 橡树 岭 国家 实验 室 (Oak Ridge National 


Laboratory, ORNL) 对 多 种 熔 盐 体系 进行 了 详细 测定 ， 为 后 续 研究 提供 了 基础 数据 上 9。 但 由 于 当时 实验 条 


点 。 研 究 人 员 通 过 实验 测量 与 相 几 计 

熔 盐 换 热 器 的 设计 离 不 3 
流 和 过 湾流 域 ， 炊 盐 在 贺 
热 过 程 的 实验 数据 还 相对 匮乏 ,主要 


EA 


件 的 限制 ， 许 多 实验 数据 的 可 靠 性 不 足 。 近 年 来 ， 随 着 熔 盐 堆 研 究 的 复兴 ， 熔 盐 热 物性 再 次 成 为 关注 的 焦 


H H 
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的 对 流 换 热 特性 与 经 
中 在 太阳 能 应 用 领域 ， 


XP HapXT ASIE 
(MSRE-PHX) (1。 该 换 热 器 是 
工 质 为 一 回路 中 的 燃料 


一 台 管 壳 式 熔 盐 - 熔 盐 换 
管 侧 ( 冷 侧 ) 工 质 为 二 回路 


的 方法 ， 获 得 了 更 加 可 靠 和 完善 的 熔 盐 物性 数据 站。 
于 对 熔 盐 传 热 与 流动 特性 的 深入 认识 。 已 有 的 实验 和 数值 模拟 结果 表明 ， 在 汕 
关联 式 预 测 吻合 较 好 ， 但 关于 熔 盐 换 热 器 内 部 复杂 传 


工 质 以 运行 温度 范围 较 低 的 HITEC JR i 73 


行 ， 没 有 出 现任 何 传 热 性 能 下 降 或 发 生 故 障 的 现象 。 


遗憾 的 是 ， 尽管 MSRE M H HREH 
里 ， 以 MSRE-PHX 为 对 象 的 深入 五 
技术 文档 、Onshape 和 Simscale 557 


熔 盐 换 热 器 性 能 提供 


了 


运行 的 熔 盐 - 熔 盐 换 热 器 的 实验 报道 最 为 详尽 的 是 MSRE 项 目 中 的 主 换 热 器 
内 器 ， 运 行 温度 范围 为 538-663°C。 其 壳 侧 ( 热 侧 》 
FP 的 冷却 盐 ， 在 超过 13000 小 时 的 实验 期 间 稳定 运 


技术 的 发 展 提供 了 重要 数据 , 在 MSRE 项 目 结束 后 的 半 个 多 世纪 
报道 。 近 期 ， 来 自 瑞 ] 
于 源 软 件 对 MSRE-PHX 进行 了 建 模 和 数值 模拟 ， 研 究 了 计算 流体 力学 
(Computational Fluid Dynamics, CFD) 技术 在 熔 盐 传 热 行为 预测 中 的 适用 性 ， 为 工业 界 采用 数值 方法 模拟 


有 益 参考 ， 但 也 指出 MSRE 历史 文献 中 的 相关 数据 存在 较 大 的 不 确定 性 。 此 外 ， 国 


的 Malcolm Aknerl' 5 基于 MSRE-PHX 的 


内 外 学 者 对 熔 盐 堆 换 热 器 的 实际 运行 特性 与 模拟 分 析 研 究 还 包括 : 中 山大 学 的 He 等 09 对 HITEC 盐 在 无 折 


流 板 管 沉 式 换 热 器 中 的 过 从 
盐 在 壳 侧 的 传 热 系数 是 Sie 


从 而 增强 了 传 热 性 能 。 
工 质 开 展 了 翅 片 管 式 众 盐 -气体 换 热 器 ， 圆 管 熔 盐 - 气 体 换 热 器 ， 以 及 管 壳 式 熔 盐 - 炊 盐 换 热 器 的 原型 试验 件 


| 流动 进行 了 实验 ，HITEC 盐 的 时 


TESTE T 


的 传 热 实验 研究 , 研究 结果 表明 管 侧 HITEC iiS m Ui 
良好 ， 层 流 区 间 实 验 数据 人 
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中 诺 数 控制 在 400-2300 内 ， 实 验 结果 表明 ， 熔 
der-Tate 关联 式 预测 值 的 3~5 倍 , 其 原因 可 能 是 壳 侧 管束 结构 形成 了 薄 的 边界 层 ， 
基 熔 盐 堆 系统 设计 和 建造 了 熔 盐 换 热 器 实验 系统 ， 以 HITEC 盐 为 


区 间 的 对 流传 热 特性 与 传统 经 验 关 系 式 吻合 
高 差 超 过 30%; M HITEC 盐 的 对 流传 热 特 性 与 导热 油 则 存在 明显 的 差别 。 陈 
玉 爽 等 (821 基 于 Ansys workbench 协同 仿真 软件 平台 , 对 比 了 300kW 的 U 型 管 式 熔 盐 换 热 器 实验 数据 和 仿 
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真 结 果 ， 并 采用 流 热 固 耦 合 方法 分 析 了 的 坎 度 场 与 应 力 场 ， 验 证 了 CFD 流体 仿真 结果 的 准确 性 。 

总 的 来 说 ， 关 于 熔 盐 堆 中 换 热 器 的 研究 仍然 较为 有 限 ， 特 别 是 在 其 传 热 与 流动 特性 的 模拟 分 析 方 面 。 
本 研究 则 在 通过 应 用 现代 计算 工具 ， 如 计算 流体 力学 软件 和 换 热 器 设计 程序 (如 HTRI Xchanger Suite), 
结合 传统 的 换 热 器 设计 方法 (如 Kern 方法 和 Bell-Delaware 方法 )， 对 MSRE-PHX 进行 详细 的 传 热 与 流动 
特性 分 析 。 本 文 将 比较 这 些 方法 的 计算 结果 与 MSRE 的 历史 运行 数据 ,验证 这 些 工具 在 高 温 炊 盐 环境 下 的 
适用 性 和 准确 性 。 


1  MSRE-PHX 设计 和 运行 参数 
11 ”结构 参数 


图 1 为 美国 橡树 岭 国 家 实验 室 10MWt 熔 盐 堆 用 的 主 换 热 器 MSRE-PHX 的 示意 图 。 其 原始 设计 于 1961 
年 完成 ， 设 备 制造 于 1963 FER, T 1964 年 进行 了 部 分 改造 ， 此 后 即 作 为 MSRE 系统 的 关键 设备 正式 运 
行 。MSRE-PHX 采用 经 典 的 管 沉 式 换 热 器 结构 ， 其 这 体 与 水 平面 呈 上 倾 3” 布 置 ， 这样 可 以 利用 重力 使 燃 
料 盐 从 壳 侧 排出 ， 同 时 也 有 利于 换 热 器 的 安装 和 维护 ， 甚 管 侧 采用 U 形 管束 ， 结 构 简 单 紧 凑 、 热 补偿 性 能 
好 。 管 束 上 加 装 了 单 己 形 折 流 板 ， 以 增加 换 热 面积 ， 提 高 换 热 效率 ， 同 时 也 减少 了 管束 的 振动 和 应 力 。 结 
构 参 数 见 表 1 。 


Inlet of Tube 
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Tube Sheet 


Inlet of Shell 


' Impingement Baffle 


OutletofTube Outlet of Shell Tube Bundle 
图 1 MSRE-PHX 的 结构 示意 图 
Fig.1 Structure diagram of MSRE-PHX 


表 1 MSRE-PHX 的 结构 参数 
Table 1 Structural parameters of MSRE-PHX 


参数 Parameter 值 Value 


折 流 板 间 距 Baffle spacing /mm 304.8 
折 流 板 厚 度 Baffle thickness / mm 6.4 

总 壳 长 Overall shell length / m 2.4384 
壳 侧 内 径 Inner diameter of shell /mm 406.4 
平均 管 长 Average tube length /m 4.2672 
管 数 Tube number 159 
管道 间距 Tube pitch / mm 19.685 
管道 外 径 Outer diameter oftube /mm 12.7 
管道 厚度 Tube thickness / mm 1.067 
有 效 换 热 面积 Effective heat transfer surface / m? 25.9 
管 板 厚度 Tube sheet thickness / mm 38.1 


12 ”设计 工 况 与 运行 参数 

MSRE-PHX 设计 工 况 P3 为 : 燃料 盐 CLiF-BeFo-ThF4-UF4-ZrF4, 70-23-1-1-5 mol?6) 以 662.78 °C 的 温度 
进入 换 热 器 后 被 冷却 到 635 *C， 体 积 流 量 为 0.075708 m3.s-1; 冷却 盐 CLiF-BeF;, 66-34 mol?) UJ 551.67 °C 
的 温度 进入 换 热 器 后 被 加 热 到 593.33 *“C， 体 积 流量 为 0.053627 m-s! . MSRE-PHX 的 运行 功率 是 由 燃料 
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盐 系 统 和 冷却 盐 系 统 的 热平衡 计算 来 确定 的 ，! 


于 在 设计 阶段 采 


传 热 性 能 不 如 预期 , 限制 了 MSRE 的 全 


的 熔 盐 物性 
功率 输出 (设计 值 为 10MW )。 表 2 是 MSRE-PHX 的 部 分 运行 参数 ， 


E BOX 


H2, MSRE-PHX 的 


表 3 是 MSRE-PHX 设计 阶段 使 用 的 物性 和 实际 物性 的 对 比 呈 ,其 中 ,当前 冷却 盐 FLiBe 的 推荐 值 由 文献 中 
中 获得 ， 将 在 后 续 计 算 中 使 用 。 


表 2 MSRE-PHX 的 运行 参数 
Table2 Operating parameters of the MSRE-PHX 


参数 Parameter 


壳 侧 Shell side 


fs 
d 
Fi 


侧 Tube side 来 源 Source 


入 口 温度 Inlet temperature ? / °C 662.78 551.67 ORNL-TM-2098P?l, 1968 
出 口 温度 Outlet temperature / °C 635 593.33 ORNL-TM-2098, 1968 
质量 流量 Mass flow * */ kgs ! 171.24 105.75 ORNL-TM-2098, 1968 
压 降 Pressure drop € / kPa 165.5 199.9 ORNL-TM-2098, 1968 
单 侧 换 热 系数 Heat transfer coefficient of each side ^/ Wm? ‘K7! 8499.9 11293.46 ORNL-TM-3039?4l, 1973 
总 换 热 系数 Overall Heat transfer coefficient */ Wm ?:K^! 3724.74 ORNL-TM-3039, 1973 
a: 设计 值 
b: 由 设计 值 
c: 测量 值 
表 3 MSRE-PHX 中 的 熔 盐 热 物性 
Table3 Thermal properties of molten salts in MSRE-PHX 
燃料 盐 冷却 盐 燃料 盐 冷却 盐 冷却 盐 
热 物性 原始 值 原始 值 修正 值 修正 值 推荐 值 
Thermal property Original values Original values Corrected values Corrected values Recommended values 
of fuel salt of coolant salt of fuel salt of coolant salt of coolant salt 
4.76 6.06 1.44 1.14 1.1 
Thermal conductivity / W-m- -K-! 
粘度 
; . 74 8.27 7.73 9.76 0.16exp(3755/ T[K]) 
Viscosity / cP 
. 2471.66 1922.22 2261.81 1971.88 2413—0.488x T[K] 
Density / kg:m-! 
比热容 
1924.64 2384.88 1981.12 2414.17 2386 
Specific capacity / J-kg- K- ! 
在 修正 了 燃 盐 物性 数据 后 ， 其 热平衡 计算 得 出 的 名 义 满 功率 为 8MWt。 但 是 ， 研 究 人 员 在 分 析 了 燃料 
盐 中 铀 和 人 同 位 素 组 成 的 变化 后 ,发 现实 际 功 率 比 SMWt 的 满 功率 要 低 799-1096. 在 当时 的 研究 中 , ORNL 
尝试 对 这 部 分 功率 偏差 做 出 解释 ， 认 为 最 可 能 的 偏差 来 源 是 冷却 盐 的 流量 测量 问题 请]。 
在 上 世纪 60 年 代 的 实验 条 件 下 , 冷却 盐 的 流量 测量 出 现 误差 是 可 以 理解 的 。 除 此 之 外 , 对 MSRE-PHX 
运行 性 能 的 观测 和 分 析 表 明 ， 传 统 的 传 热 关系 式 也 适用 于 熔 盐 传 热 。 在 超过 14000 小 时 的 运行 时 间 内 ， 在 
MSRE 的 运行 条 件 下 ， 熔 盐 没 有 对 换 热 器 的 结构 材料 CHastelloy-N 合金 ) 产生 明显 的 腐蚀 或 者 结 垢 现象 ， 


因为 实际 观测 表明 MSRE-PHX 的 传 热能 力 术 


因 腐 蚀 、 结 


构 设 计 和 材料 选择 都 是 合 


MSRE-PHX 的 结 


E 和 有 效 的 。 在 设计 


MSRE-PHX 实际 传 热 性 能 不 如 预期 的 主要 原因 ， 而 冷却 盐 的 流量 
低 于 名 义 功 率 的 最 大 可 能 原因 。 

2 ”研究 方法 

2.1 理论 计算 方法 

管 壳 式 换 热 器 的 理论 计算 可 以 分 为 管 侧 和 壳 侧 两 方面 。 本 文中 ， 
有 折 流 板 的 壳 侧 采用 两 种 计算 方法 : Kern 法 和 Bell-Delaware 法 71 
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后 、 旁 路 或 流动 限制 而 下 降 。 这 说 明 


阶段 使 月 


DH 


日 有 较 大 1 


则 计算 采 月 


遍 差 的 熔 盐 热 物 性 数据 是 
测量 误差 是 MSRE-PHX 的 实际 传 热 功 率 


中 Kern 


N 


H Gnielinski 方法 PC9， 带 
法 诞生 于 1950 年 ， 是 
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MSRE-PHX 在 设计 阶段 采用 的 方法 ，Bell-Delaware 法 诞生 于 1963 年 ， 是 目前 公开 发 表 的 各 种 方法 中 最 先 
进 的 一 种 。 
2.1.1 Mit- Kern 法 
Kern 法 基于 采用 标准 公差 制造 的 商用 管 壳 式 换 热 器 的 大 量 实验 研究 ， 是 一 种 经 典 的 换 热 器 设计 方法 ， 
它 的 主要 优点 是 简单 易 用 ， 缺 点 是 计算 壳 侧 传 热 系数 和 压 降 时 忽略 了 壳 侧 的 漆 露 流 影响 。 
由 于 壳 侧 的 流动 截面 积 沿 壳 体 长 度 变 化 ，Ker 法 中 的 壳 侧 流速 w 根据 壳 侧 错 流 面积 4 计算 。 壳 侧 的 


LES 


等 效 直径 D, 则 随 管束 排列 方式 不 同 ，MSRE-PHX 中 的 管束 排列 为 等 边 三 角形 ， 因 此 : 


= 于 
2 
4| 至 0.87P -0.51% 
2 4 


x 元 (1) 
A—— 
2 
壳 侧 错 流 面 积 4 : 
A, s (È E d, )D,l, ; Q) 
' P 
" EN EX n. TR G. : c cas G.D Dop. ; 
因此 可 以 分 别 计算 出 壳 全 质量 流速 G, = 站 和 流速 上 = 二， 以 及 雷诺 数 Re= o Ese Kem 法 
S S S Hs 
中 的 努 塞 尔 数 Nu, 表达 形式 为 
hD 0.14 
Nu, = a 6) 
对 于 壳 侧 压 降 ，Kem 法 同样 给 出 基于 壳 侧 雷诺 数 的 无 量 纲 摩擦 系数 j; ， 因 此 : 
"EL L DV; H 、-014 
AP =$ S S . 4 
: (2E 2 PE S 


24.2 ” 壳 侧 计算 -Bell-Delaware 法 

Bell-Delaware 法 中 考虑 了 影响 壳 侧 流体 分 布 的 诸多 因素 ， 更 为 精确 。 其 设计 方法 以 Tinker 提出 的 流动 
模型 为 基础 ， 将 壳 侧 流体 分 为 5 股 流 路 ， 如 图 2 所 示 。 

这 5 股 流 路 分 别 为 : 
流 路 A: 折 流 板 管 孔 与 管子 之 间 的 泄露 流 ; 
流 路 B: 错 流 ， 与 不 考虑 泄露 流 路 的 理想 流 路 方向 相同 ， 是 影响 壳 侧 传 热 和 流动 阻力 的 主要 流 路 ; 
流 路 C: 管束 外 围 和 壳 内 壁 之 间 的 旁 流 ; 
WERE: PTAA SE A EE Z E E ER Mik; 
流 路 F: 管 程 分 程 隔 板 中 间 位 置 的 旁 流 。 
在 设计 管 壳 式 换 热 器 时 ， 为 了 获得 更 好 的 传 热 效果 ， 流 路 了 B 的 占 比 应 越 大 越 好 。 


tb 
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Hot Fluid Cold Fluid PPT " 
Inlet Inlet i BaMe 1 Shell Rear-End 
A 


图 2 管 壳 式 换 热 器 壳 侧 流 路 示意 图 [28] 
Fig.2 Flow diagram of shell side in a STHE 


l. seu S EIRA ERA EAP LASAO EAS P: 
HURAR S, 可 以 表示 为 : 


S, =l, |o. -D,) ea p ZI 


n 


其 中 ， 书 是 垂直 于 流动 方向 的 管 间距 。 
穿 过 管束 最 外 层 管 子 中 心 的 假想 圆 的 直径 D。 : 


Du -D, -d, 


折 流 板 切 口 与 管束 直径 相交 对 应 的 中 心 张 角 0,: 


ctl 


0, = zi 


其 中 B. 是 缺口 高 度 与 党 体内 径 的 比值 。 


(5) 


(6) 


(7) 


HARER KETAT ETAM: F, =- (0, -sing) ;两 折 流 板 缺 口 区 之 间 的 错 流 区 中 的 管 
pA 


子 数 与 管子 总 数 的 比值 FL-1-2F, - 14 (sin6, -0 )。 
T 


v A MN S y > Av A ME D 1-2B 
两 折 流 板 缺 口 区 之 间 的 钳 流 区 内 错 流 方 向 上 的 管子 行 数 ，N, o 0C 202 
P 
ETTI N, =D., 
P 


; 缺口 


区 内 暴露 于 错 流 的 有 


上 式 中 ， 由 于 缺口 区 的 流动 部 分 平行 于 管子 轴线 ， 部 分 垂直 于 管子 轴线 ， 因 此 取 系 数 0.8。 忆 是 错 流 


方向 上 的 管 间 距 。 
对 于 流 路 A， 即 折 流 板 管 孔 与 管子 之 间 的 泄露 流 的 面积 So 可 由 下 式 计 算 ; 


S, = EE +25 -42) Nar) 
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(8) 
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56, 是 折 流 板 孔 与 管 外 径 之 间 的 直径 间隙 ， 未 知 时 可 按 TEMA 标准 取 0.8 或 者 0.4。 
对 于 流 路 忆 ， 即 折 流 板 与 这 内 壁 之 间 的 泄露 流 ， 其 泄露 面积 可 以 表示 为 


S, — Dó,, Gr E 0.50,.) (9) 


其 中 ，0 =2cos"(1-2B.) ; 6, 是 折 流 板 直 径 与 党 内 径 之 间 的 直径 间隙 ， 可 按 此 式 计算 : 
ô, =0.8+0.002D 。 
对 于 流 路 C， 即 管束 与 壳 内 壁 之 间 的 泄露 流 ， 其 截面 积 为 ; 


S, -1, (D, -D,) (10) 


对 于 流 路 B， 即 错 流 ， mn 区 的 净 错 流 截面 积 为 总 错 流 截 面积 减 去 管子 所 占 的 截面 积 之 值 : 
S. l5 (0, — sinQ,.) LN Fund? (11) 


在 Bell-Delaware 方法 中 ,为 了 更 接近 实际 流动 ， 首 先 要 计算 壳 侧 流动 垂直 于 管束 的 理想 错 流 条 件 下 的 
传 热 系数 和 摩擦 系数 ， 然 后 将 理想 情况 下 所 得 结果 乘 以 根据 流 路 校正 因子 ， 从 而 得 到 最 终结 果 。 
过 侧 传 热 系 数 为 ; 


h, = hea ffijyjs (12) 
XB. A, 是 理想 管束 纯 错 流 情 况 下 的 壳 侧 传 热 系 数 : 


À 2/3 0.14 

H m, s s 

ho = J idealC p,s E J: u | E | (13) 
m pss SW 


1.33 Y 
其 中 , Tea = alza] Re” |» Qq— ^ Q0 O55 05 和 a, 取 自 相 关 经 验 表 格 。 
e. ' 


Bell-Delaware 方法 中 的 特征 长 度 为 管 外 径 4, 。 壳 侧 质量 流速 C, =m, /9, ， 壳 侧 流速 ww = Gp, , se 


_ pv.d, d, 
Hs 


2， 流 路 传 热 校正 因子 的 计算 : 
折 流 板 切 口 和 跨 距 校正 因子 j.: 


雷诺 若 数 Re, 


j, 20.554 0.72F, (14) 


折 流 板 与 管子 、 壳 内 壁 之 间 泄 露 流 的 校正 因子 j : 
=0.44(1—x.)+[1—0.44( -r ]exp(72.2r,,,) (15) 


r= So r = S sb Ẹ S, tb 
s S 十 S ? Im S 2 
sb tb 


m 


管束 流 路 C 和 流 路 下 的 旁 路 校正 因子 广 : 


S 
dej]. cus 


m 


N 
4 
X 


(16) 
1 r, 20.5 


KH, n-NI/N,  Re<100 时 ，C, =1.35 ， 当 Re>100 时 ，C, =1.25 ，N, 为 密封 挡 板 对 数 ， N. 为 


两 块 折 流 板 之 问 的 错 流 管 排 数 。 
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进出 口传 热 校正 因子 


ro (N, 7 1) F (L; /le i» eu (L, / Lo ya 
! (N, - D (1, ZELUM his) 


其 中 ， 当 Re<100 Hf, 5 =0.6， 当 Re>100 Bf, n =1/3。 
层 流下 逆向 传 热 温 差 校正 因子 记 ， 仅 在 Re<100 时 才 做 考虑 ， 本 研究 中 并 不 涉及 ， 因 此 取 值 为 1， 不 
HAR. 

3， 压 降 计 算 

Bell-Delaware 方法 中 将 壳 侧 压 降 分 为 三 部 分 ， 分 别 是 端 部 区 〈 进 出 口 区 ) 错 流 管束 压 降 AP 、 非 端 部 
区 错 流 管束 压 降 AP 和 人 缺口 区 管束 压 降 AP, ， 即 : 


(17) 


AP, - AP. + AP, AF, (18) 


非 端 部 区 错 流 压 降 AP 的 计 


AR = (N, E DAE, RR, (19) 


式 中 ，AP, 是 理想 情况 下 的 错 流 压 降 : 


f N G? 0.14 
AP al = 2 ideal c —s [es (20) 
8P LH, 
iay b 
Jt. fau b AE) Reh, bm oss So bo LR DG n 


Jb 5h. g 是 单位 转换 因子 ， 值 为 1; RR 是 泄露 流 路 校正 因子 ， R -exp[-1.330 7,7] ， 


p=0.8-0.15(1+r,); RR, 是 旁 流 校正 因子 ， 当 ,<0.5 时 ， R, =exp| -C, (5,58, (1-427... 否则 R，=1。 


对 于 Re>100，C，=3.7。 


对 于 缺口 区 压 降 ， 当 Re>100 时 ， 理 想 情 况 下 的 压 降 : 


240.6N n? 
Lo a QU) 
2g.05,8,, 
因此 ， 实 际 缺 口 区 压 降 : 
AP, 一 N,AP, idea Py (22) 
对 于 端口 区 压 降 : 
N. 
AP, = 2AP ral í + 2d R,R, (23) 


Ep. R= os 4. 75,) O D UNT Jem L 和, 分别 是 折 流 板 进出 口 区域 的 距离 。 


2.2 HTRI Xchanger Suite 


HTRI Xchanger Suite 是 由 美国 传 热 研究 协会 (Heat Transfer Research Institute， HTRI) 开发 的 专业 换 
热 器 设计 、 校 核 和 模拟 软件 。 该 软件 基于 半 个 多 世纪 以 来 收集 的 工业 传 热 设备 试验 数据 ， 能 够 灵活 且 严 格 
的 设计 换 热 器 的 几何 结构 ， 从 而 利用 HTRI 独 有 的 经 验 公式 对 换 热 器 的 性 能 进行 精确 评估 [231。 
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换 热 器 的 理论 计算 根据 所 需 工 况 ， 使 用 通用 的 传 热 方程 和 物性 数据 来 估算 换 热 器 的 性 能 ， 并 不 考虑 换 
热 器 内 部 的 具体 情况 。 而 HTRI 作为 专门 的 换 热 器 设计 计算 软件 ,考虑 了 具体 的 换 热 器 设计 参数 ， 如 管束 
布局 、 管 束 尺 寸 、 流 体 流动 等 ， 并 将 管 程 和 壳 程 分 成 多 个 区 间 ， 并 对 每 个 区 间 进 行 独立 的 传 热 和 流体 力学 
迭代 计算 ， 因 此 可 以 更 准确 地 估算 换 热 器 的 性 能 ， 并 考虑 到 换 热 器 内 部 的 具体 情况 。 


23 CFD 仿真 方法 


为 了 能 够 更 好 的 了 解 MSRE-PHX 内 的 流动 和 传 热机 理 ， 本 研究 中 使 用 了 商业 CFD 软件 Ansys Fluent 
对 MSRE-PHX 进行 了 全 尺寸 建 模 和 分 析 ， 对 于 该 换 热 器 中 的 流速 、 压 力 和 温度 场 的 具体 分 布 和 对 换 热 器 
性 能 产生 的 影响 进行 了 详细 的 探讨 。 此 外 ， 本 研究 中 也 计算 了 传 热 系数 和 压 降 等 全 局 参数 ， 并 与 基于 经 验 
关系 式 的 理论 计算 和 专业 换 热 器 设计 软件 HTRI 进行 比较 分 析 , 探索 了 对 大 型 换 热 器 进行 全 尺寸 CFD 建 模 
和 分 析 的 可 能 性 和 局 限 性 。 

2.3.1 ”计算 模型 及 物性 参数 

计算 模型 如 图 1 所 示 ， 计 算 所 用 物性 参数 如 表 3 所 示 。 

2.3.2 ”边界 条 件 

数值 模拟 中 的 边界 条 件 设置 如 下 : 

入 口 边界 : 两 入 口 均 为 质量 流量 入 口 。 

出 口 边界 : 压力 出 口 。 

壁面 边界 : 壳 侧 外 壁面 为 无 滑 移 的 绝热 壁面 。 

2.3.3 ”网 格 划 分 

MSRE-PHX 的 计算 域 包括 燃料 盐 流 体 域 、 冷 却 盐 流体 域 和 其 余 金 属 结构 组 成 的 固体 域 。 

由 于 MSRE-PHX 尺寸 大 〈 约 为 2.5m*0.5m*0.5m)、 管 束 多 (159 JR U 形 管 )， 在 经 多 次 试验 后 ， 本 研 
究 中 采用 商业 软件 Fluent Meshing 对 模型 进行 网 格 划 分 。Fluent Meshing 的 特点 是 具有 网 格 划分 速度 高 ， 网 
格 处 理 灵 活 度 高 ， 并 且 能 够 生成 基于 “马赛 克 ” 网 格 技术 的 Poly-Hexcore 体 网 格 ， 尽 可 能 提升 体 网 格 中 六 面 
体 的 数量 ， 从 而 提升 求解 效率 和 精度 B91。 

为 尽 可 能 准确 模拟 熔 盐 在 换 热 器 内 的 流动 ， 在 过 体内 壁 、 管 束 等 流体 域 与 固体 域 交界 处 设置 边界 层 网 
格 。 对 于 边界 层 网 格 ，Fluent Meshing 同样 可 以 运用 “马赛 克 ” 技 术 ， 实 现 壁 面 附 近 的 层 状 Poly 网 格 、 过 渡 
区 域 的 纯 Poly 网 格 和 核心 区 域 的 六 面体 网 格 过渡 , 进一步 提 升 网 格 整体 质量 ， 且 有 效 降低 网 格 总 体 数量 与 
求解 时 间 。 网 格 示意 图 如 图 3 所 示 。 


(b) (c) 


图 3 MSRE-PHX 的 网 格 划分 示意 图 
Fig.3 Meshing of MSRE-PHX 
2..4 ”网 格 无 关 性 分 析 和 求解 控制 设置 
网 格 划 分 数量 分 别 为 38957352、55713541 和 74893211 时 ， 以 迭代 残 差 和 质量 流量 、 温 度 平衡 为 收敛 
条 件 ， 计 算 迭 代 约 5000 次 后 达到 收敛 ， 壳 侧 出 口 温度 分 别 为 637.5、639.61 和 639.4*%C， 综 合 考虑 计算 资 
源 和 时 间 ， 选 择 网 格 数量 为 55713541 时 ， 和 迭代 10000 次 的 计算 结果 进行 分 析 。 
在 计算 时 选择 COUPLED 算法 来 求解 ， 该 算法 可 用 更 少 的 迭代 次 数 收 敛 。 计 算 收 敛 的 判断 准则 是 各 项 
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物理 量 的 迭代 残 差 ， 为 保 说 
口 的 质量 流量 和 温度 ; 


结果 和 讨论 


3.1 理论 计算 校 核 


3 


值 为 3724.74W.m-2.K-1，Bell-Delaware 关系 式 计 


值 ，Kern 关系 式 
数 ，Gnielinski 计 


论 换 热 系 数 偏 高 。 由 于 Gnielinski 关系 式 中 
强 。 在 壳 侧 ， 两 种 关系 式 都 是 对 壳 侧 流体 进行 了 理想 化 假设 ， 导 致 换 热 系数 的 理论 
可 结构 带 来 的 流 路 变化 ， 偏 于 保守 ， 因 


系 式 中 并 未 考虑 壳 侧 几 


结果 


B BJ] JA IB ALS 


表 4 是 MSRE-PHX 的 理论 校 核 结果 。 在 传 热 系数 方面 ，MSRE-PHX 实验 所 测 


得 到 的 总 换 热 系数 为 3506.3W:m-?:K-!， 
得 出 的 管 侧 换 热 系数 与 实验 测量 
式 计 算 的 壳 侧 换 热 系 数 与 实验 测量 值 相 比 则 分 别 相差 4.7% 和 62%， 即 


到 的 


得 


直 相 差 在 -14% 以 内 


不 考虑 弯 管 内 的 换 热 ， 


式 考虑 的 相对 全 面 ， 但 是 在 MSRE-PHX 的 实验 中 ， 实 际 
里 想 假设 差距 很 远 ， 以 及 本 文中 管 侧 冷却 盐 
的 不 同 ， 也 导致 了 Bell-Delaware 关系 式 的 计算 得 


置 的 隔 热 层 与 


在 压 降 方面 ， 壳 侧 压 降 的 实验 值 为 165.5 kPa. EN 


实验 室 测 量 流量 的 方式 是 使 用 基于 压 差 的 文 丘 里 管 流 量 计 
管 侧 流 量 减 小 导致 管 侧 压 降 减 小 ， 


当 仅 作为 参考 。 在 表 4 


中 ， 


与 实验 值 的 偏差 在 可 以 接受 的 范 


Bell-Delaware 关系 式 计 


EA. EFEN 
的 壳 侧 压 降 明显 
表 4 


的 结构 参数 如 U 形 管 


x 
到 


IHE, 
Wm (~250kPa)。 


1 


上 BERIN 


Kern 关系 式 计算 的 壳 侧 月 


完全 收 公 ， 对 和 迭代 残 差 标准 统一 设置 为 小 于 102。 除 了 监视 迁 代 残 差 ， 对 j 
行 同步 监视 ， 以 加 


2H 
F 


的 总 换 热 系数 的 平均 
总 换 热 系数 为 4087.3W-m-?.K-!, 高 于 实验 测量 
氏 于 实验 测量 值 。 比 较 换 热 器 两 侧 各 自 的 换 热 系 
>, Kern 关系 式 和 Bell-Delaware 关系 


ams des 


WEE 


烈 


IW? 


V3 iru AERE Ji 


In] 


管 侧 的 理论 换 热 系数 偏 低 ,， 壳 侧 的 理 


实际 传 热 增 


实际 值 。Kern 关 


此 计算 的 换 热 系数 较 低 ，Bell-Delaware 关系 


三 | 
FH. 


fs 


zx] 


MSRE-PHX 的 理论 校 核 结 果 


Table4 Theoretical verification results of MSRE-PHX 


的 排 布 方式 、MSRE-PHX 设 
的 物性 值 也 与 当时 实验 计算 换 热 系数 时 使 用 
的 换 热 系数 高 于 实际 值 。 
I 压 降 的 实验 值 为 199.9 kPa。 需 要 说 明 的 是 ,橡树 岭 
测量 有 误 ， 那 么 压 降 的 实验 测量 值 应 
党 侧 压 降 的 计算 值 分 别 为 ~211 和 ~177 kPa, 
E 降 明显 偏 低 (~140kPa)， 


管 侧 Tube side 壳 侧 Shell side 管 侧 Tube side FEM Shell side 


参数 Parameter 


Gnielinski Kern Gnielinski Bell-Delaware 

出 口 温 度 Outlet temperature / °C 577.66 643.45 580.91 641.03 
压 降 Pressure drop /kPa 211.82 140.09 211.5 250.26 
单 侧 换 热 系 数 

AE . usi 9742.9 8901.9 9808.1 13771.4 
Heat transfer coefficient of each side / W:m ^K 
总 换 热 系数 Heat transfer coefficient of each side / W:m? K! 35063 4087.3 
总 换 热 量 Total heat transfer rate / MW 6.55776 7.3774 


32 HTRI 模拟 计算 结果 


图 4 (a) - Ce) 展示 了 MSRE-PHX 在 HTRI 中 的 分 段 计算 结 果 。 图 4 (a) fH 


体 和 


职 温度 治 长 度 方 向 的 逐步 变化 情况 。 


Z] 


压 降 变化 趋势 总 体 平稳 


给 出 了 两 侧 流体 在 流动 过 程 中 所 传递 
过 程 中 的 Re 数 ， 可 以 看 出 除 
两 侧 换 热 系数 和 总 换 热 系 数 沿 长 度 方 
16536W:m-2.K-1 


系 数 为 


， 但 在 壳 侧 入 


了 壳 侧 的 入 


113kPa， 
fi 


E 


E 


, 


总 换 A 
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4 CO 中 为 管 程 和 壳 程 流体 在 每 个 
处 的 压 降 达 到 
的 累计 热 内 总 量 为 8.1222MW，。 
口 段 和 出 口 段 ， 主 要 换 热 长 度 


可 的 变化 ， 最 终 输 


又 间 的 压 降 
区 间 的 压 降 均 低 于 40kPa。 


Z] 


A 


ZIN 


其 


上 了 管 程 和 壳 程 流体 的 
直 ， 可 以 看 出 两 侧 


图 4 Ce) 中 


4 d) 中 为 两 侧 流体 在 流动 


区 间 中 Re 值 变化 不 大 。 


系 数 为 


Z] 


4 Ce) 为 


的 管 侧 换 热 系数 为 8012.1 W-m-?-K-',. 5l f 2A 
3731.1W.m-2.K-1 


作者 


: 中 文 题目 
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的 模拟 结果 


Fig.4 Simulation results of HTRI 


33 CFD 仿真 结果 


3.3.1 


图 5 为 MSRE-PHX 的 MSRE-PHX 的 温度 分 布 云 


fiL TE 


4? 


燃料 盐 


温度 场 分 析 


Z] 


nm 


EMD 和 温度 分 布 流 线 图 ( 管 侧 ), 使 月 


2400 


真 所 得 出 的 壳 侧 出 口 温度 为 639.61°C, 管 侧 出 口 温 
IERE U 形 管 流动 过 程 中 不 断 吸收 热量 ,温度 从 入 口 至 


度 为 582.85°C。 


Ma app 


背 流 区 ， 


R. 无 法 与 只 


管束 ， 这 种 冲刷 增强 了 


E 出 口 逐 步 升 高 。 而 在 党 侧 ， 
出 流 程度 ， 从 而 增强 


Hd 
中 


zd 


N 


温度 都 相对 较 低 。 其 原因 是 折 流 板 背 流 


又 存在 


它 区 域 的 燃料 盐 进 行 有 效 的 热量 交换 。 


行 更 多 的 热 


TN 


量 交换 ， 从 而 导致 其 温度 相对 低 于 其 它 区 1 


B--a 


回流 现象 ， 在 这 个 区 域 ， 


的 燃料 盐 


管 侧 冷却 盐 的 温度 分 布 趋势 明显 ， 
ET UR EGET DL 
一 块 折 流 板 
燃料 盐 的 流动 受到 阻 
相反 ， 这 部 分 燃料 盐 会 与 管束 以 及 管束 内 的 冷却 盐 进 
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图 5 MSRE-PHX 的 温度 分 布 云图 与 流 线 图 
Fig.5 The cloud map of temperature distribution and flow paths colored with temperature of MSRE-PHX 

3.3.2 ”压力 场 分 析 

图 6 为 MSRE-PHX 的 压力 分 布 云图 〈 壳 侧 ) 和 压力 分 布 流 线 图 ( 管 侧 )。 由 于 折 流 板 的 存在 ， 燃 料 盐 
在 冲击 折 流 板 后 流向 发 生 显著 变化 ， 形 成 了 较 大 的 速度 梯度 。 这 是 由 于 折 流 板 的 阻挡 和 扰 流 作用 改变 了 流 
体 的 流动 方向 和 速度 ， 从 而 引起 了 速度 的 突变 和 梯度 的 产生 。 此 外 ， 当 燃料 盐 通 过 折 流 板 时 ， 其 流通 截面 
积 会 发 生 显著 变化 , 这 会 导致 燃料 盐 的 流速 发 生 相 应 的 变化 。 这 种 流速 的 突变 也 会 对 流体 的 压力 产生 影响 ， 
从 而 产生 较 大 的 局 部 压 损 。 在 折 流 板 的 背面 , 由 于 流体 的 阻挡 和 扰动 , 会 形成 一 定 的 回流 区 。 在 这 个 区 域 ， 
流体 的 流动 方向 和 速度 都 会 发 生变 化 , 这 也 会 形成 局 部 压 损 。 折 流 板 的 向 流 和 背 流 处 的 压 差 对 折 流 板 的 结 
构 强 度 提 出 了 一 定 的 要 求 。 如 果 折 流 板 的 结构 强度 不 足 ， 可 能 无 法 承受 这 种 压 差 所 产生 的 力 ， 从 而 导致 折 
流 板 的 损坏 。 因 此 ， 在 设计 折 流 板 时 需要 充分 考虑 其 结构 强度 ， 以 确保 其 能 够 在 换 热 器 的 运行 过 程 中 保持 
稳定 和 安全 。 
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图 6 MSRE-PHX 的 压力 分 布 云图 与 流 线 图 
Fig.6 The cloud map of pressure distribution and Flow paths colored with pressure of MSRE-PHX 

3.3.3 ”速度 场 分 析 

图 7 为 MSRE-PHX 的 速度 分 布 云图 〈 壳 侧 ) 和 速度 分 布 流 线 图 ( 管 侧 )。 从 图 中 可 以 看 出 ， 单 弓形 折 
流 板 的 存在 使 MSRE-PHX 的 壳 侧 燃料 盐 呈 现 明 显 的 周期 性 “Z” 形 流动 。 在 折 流 板 的 圆 缺 处 ， 燃 料 盐 的 流 
向 与 管束 一 致 ， 即 为 顺 流 ， 此 时 速度 梯度 与 温度 梯度 方向 垂直 ， 根 据 场 协同 理论 ， 此 时 换 热 效果 不 佳 ; 在 
折 流 板 的 向 流 区 ， 燃 料 盐 横 向 冲刷 管束 ， 即 为 错 流 区 ， 此 时 速度 梯度 与 温度 梯度 方向 近乎 一 致 ， 此 时 换 热 
效果 最 好 。 
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图 7 MSRE-PHX 的 速度 分 布 云图 


Fig.7 The cloud map of velocity distribution and flow paths colored with velocity of MSRE-PHX 


图 8 为 MSRE-PHX 的 速度 矢量 图 ( 壳 侧 )。 可 以 看 出 ,燃料 盐 流 过 折 流 板 后 ， 背 流 区 的 速度 明显 降低 ， 
如 图 8 Ca) 所 示 ， 折 流 板 背 部 可 见 明 显 回流 区 。 此 处 燃料 盐 近 乎 停滞 ， 形 


对 该 区 域 的 速度 矢量 放大 观察 ， 


成 所 谓 的 流动 “ 死 区 ”， 可 能 会 产生 较 严 习 


的 结 折 沉 积 问题 ， 这 也 与 温度 和 压力 场 的 分 布 情况 相 呼应 。 折 流 


板 的 小 缺口 形成 了 漏 流 区 ， 该 区 域 的 存在 使 得 该 处 的 部 分 流体 无 需 绕 过 折 流 板 就 能 直接 到 达 折 流 板 的 背 流 
区 ， 增 加 了 顺 流 区 域 ， 减 小 横流 引起 的 流 致 振动 几率 ， 降 低 了 壳 侧 阻力 ， 减 组 了 该 处 的 结 垢 速率 。 总 的 来 


说 ， 圆 缺 处 的 流速 高 于 错 流 区 的 流速 ， 而 背 流 
征 ， 这 种 现象 会 导致 较 大 的 压 损 
出 现在 燃料 盐 撞 击 防 冲 板 处 。 在 入 口 区 由 于 防 冲 板 的 作用 


区 的 流速 最 低 ， 这 是 壳 侧 装 有 单 号 形 折 流 板 时 的 典型 流动 特 


产生 ， 同 时 也 会 降低 换 热 效率 。 此 外 ， 由 图 8 Co 可 知 ， 壳 侧 的 最 大 流速 


燃料 盐 被 分 为 两 股 ， 并 产生 小 的 旋涡 ， 流 速 也 


~ 


因此 增加 。 燃 料 盐 流 过 防 冲 板 后 ， 流 速 明显 降低 ， 因 此 对 入 口 处 的 管束 不 会 产生 过 大 的 冲击 ， 从 而 降低 了 


管束 振动 现象 的 产生 。 在 燃料 盐 出 口 区 , 为 了 降低 燃料 盐 的 高 温 对 管 板 的 影响 , 设置 了 一 个 燃料 盐 滞 留 区 ， 


如 图 8〈c)。 燃 料 盐 经 一 狭窄 入 口 进入 该 滞留 区 后 ， 几 乎 停止 流动 ， 结 合 图 5 的 温度 分 布 ， 可 以 看 出 灌 留 


的 结构 完整 性 。 


区 的 燃料 盐 温度 明显 低 于 主流 区 的 燃料 盐 温 度 ， 


因此 管 板 上 不 会 产生 大 的 温度 梯度 变化 ， 从 而 保护 了 管 板 
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图 8 MSRE-PHX 的 速度 矢量 图 
Fig. 8 Velocity vector diagram of MSRE-PHX 
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3.4 结果 对 比分 析 


K 5 为 MSRE-PHX 实验 与 理论 计算 、HTRI 软件 和 仿真 分 析 对 比 的 汇总 表 。 从 表 5 可 以 看 出 ， 无 论 是 
总 换 热 系数 U 还 是 换 热 量 Q,， 两 种 理论 计算 方法 和 两 种 仿真 模拟 的 结果 与 实验 结果 的 差异 均 在 可 接受 范围 
内 。 其 中 Kern 方法 与 MSRE 实验 所 得 的 换 热量 相差 最 大 ， 为 15%，HTRI 软件 模拟 与 MSRE 实验 所 得 的 
总 换 热 系数 相差 最 小 , 为 0.16%。 由 此 表明 理论 计算 和 仿真 分 析 对 MSRE-PHX 性 能 预测 均 有 相当 的 适用 性 
和 有 效 性 。 

首先 ， 在 传 热 方面 ， 在 换 热 器 壳 侧 出 口 温 度 的 预测 上 ， HTRI 软件 和 CFD 仿真 的 计算 结果 与 MSRE 
实验 结果 的 差距 非常 小 ， 误 差分 别 为 0.6% 和 0.72%， 说 明 对 于 复杂 的 壳 侧 流动 ， 这 两 种 方法 的 预测 相当 准 
确 。 而 对 于 管 侧 的 冷却 盐 出 口 温 度 ，HTRI 软件 和 CFD 仿真 计算 值 与 MSRE 实验 结果 相 比 稍 低 ， 误 差分 别 
为 -1.6% 和 -1.7%。 

其 次 ， 在 压 降 方面 ， 对 于 壳 侧 压 降 和 管 侧 压 降 ， 理 论 计算 和 仿真 分 析 的 结果 均 比 实验 结果 高 。 
MSRE-PHX 的 报告 中 提 到 由 于 壳 侧 入 口 和 出 口 处 存在 较 大 的 压 损 ， 实 际 压 降 可 能 会 达到 实验 值 的 两 倍 ， 这 
种 现象 与 HTRI 和 CFD 仿真 的 结果 一 致 ， 因 为 在 HTRI 中 壳 侧 入 口 的 压 降 达到 113kPa， 在 CFD 仿真 中 壳 
侧 入 口 的 压 降 达 到 349.03KPa。 此 外 ， 当 时 实验 条 件 的 限制 〈 如 熔 盐 对 测量 仪器 的 腐蚀 ) 也 可 能 导致 所 测 数 
据 的 偏差 ， 导 致 实验 结果 与 理论 或 仿真 分 析 的 预测 结果 存在 差异 。 

最 后 ， 分 别 使 用 理论 计算 、HTRI 软件 和 CFD 仿真 三 类 四 种 方法 对 MSRE-PHX 的 性 能 进行 分 析 并 汇 
总 后 ， 可 以 发 现 ， HTRI 软件 模拟 和 CFD 仿真 的 结果 与 实验 数据 最 为 接近 ， 而 Kern 方法 作为 一 种 保守 的 
理论 计算 方法 , 与 实验 结果 最 为 偏离 。 这 说 明 Kern 方法 在 换 热 器 计算 中 存在 的 一 些 假设 或 简化 会 使 结果 与 
实际 运行 性 能 产生 较 大 偏差 .而 HTRI 软件 和 CFD 方法 可 以 更 好 地 模拟 复杂 的 流动 和 传 热 过 程 。 作 为 一 种 
改进 的 理论 计算 方法 ，Bell-Delaware 方法 与 同样 基于 流 路 分 析 法 的 HTRI 软件 计算 结果 较为 接近 ， 但 是 计 
算 过 程 十 分 繁琐 。HTRI 软件 和 CFD 仿真 分 析 相 对 更 为 精确 ， 却 也 存在 差异 。HTRI 软件 相 较 于 CFD 仿 
真 更 为 简便 ， 尤 其 是 在 换 热 器 结构 是 商业 上 成 熟 的 型 式 时 ，HTRI 软件 非常 适合 用 作 工 程 设计 参考 。 而 要 
深入 研究 换 热 器 内 部 的 温度 分 布 和 流动 特征 , 特别 在 采用 新 兴 传 热 工 质 、 新 型 流 场 结构 设计 新 式 换 热 器 时 ， 
需要 使 用 CFD 仿真 方法 。CFD 方法 需要 更 多 的 计算 资源 ， 但 能 够 给 出 更 加 精确 详细 和 直观 的 信息 。 

表 5 MSRE-PHX 实验 与 理论 计算 、 仿 真 分 析 对 比 


Table5 The comparison of experiment results with theoretical calculation and simulation analysis of MSRE-PHX 


参数 偏差 偏差 偏差 偏差 
MSRE Kern i Bell . HTRI j CFD : 
Parameter Discrepancy Discrepancy Discrepancy Discrepancy 


壳 侧 出 口 温 度 

. 635 - -- - - 638.85 - 639.61 - 
Outlet temperature of shell side / °C 

管 侧 出 口 温度 

Outlet temperature of tube side / °C 


壳 侧 压 降 


Fé 
EL 


593.33 - - - - 583.84 - 582.8  - 


i 165.48 192.9 16.56% 222.4 34.4% 372.26 124.96% 349.03 110.92% 
Pressure drop of shell side / kPa 


管 侧 压 降 


EY 1A) Hh. pE 


. 199.96 310.27 55.1796 310.27 55.1796 298.13 49.09% 226.54 13.29% 
Pressure drop of shell side / kPa 


总 换 热 系数 
3725 344423 -7.54% 3990.09 7.1296 3731.1 0.16% 3510.15 -5.77% 
Heat transfer coefficient of each side / W:m ?: K^! 
总 传 热量 
8 6.8 1596 7.89 -1.38% 8.12 1.5% 7.86 -1.75% 


Total heat transfer rate / MW 


本 文 以 美国 橡树 岭 国家 实验 室 10MWt 熔 盐 扒 用 的 主 换 热 器 为 研究 对 象 ， 基 于 该 换 热 器 的 设计 参数 ， 
应 用 管 壳 式 换 热 器 的 理论 计算 方法 (包括 Kern 方法 和 Bell-Delaware 方法 )、 换 热 器 设计 软件 HTRI Xchanger 
Suite 和 CFD 仿真 方法 对 换 热 器 的 关键 性 能 指标 (如 换 热 系数 、 压 降 和 换 热 功率 ) 分 别 进行 了 计算 分 析 ， 
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得 到 以 下 结论 : 

1. 首先 ， 对 两 种 理论 计算 、HTRI 软 件 及 CFD 仿 真 模拟 等 三 种 换 热 器 的 设计 方法 进行 了 详细 的 分 析 ， 三 种 
设计 方法 的 结果 与 实验 乡 果 的 差异 均 在 可 接受 范围 内 ， 其 中 Kern 方 法 与 MSRE 实 验 所 得 的 换 热 量 相差 最 
大 ， 为 15%，HTRI 软 件 模 拟 与 MSRE 实 验 所 得 的 总 换 热 系数 相差 最 小 ， 为 0.16%， 表 明理 论 计算 、HTRI 

软件 及 CFD 仿 真 分 析 对 MSRE-PHX 性 能 预测 均 有 相当 的 适用 性 和 有 效 性 

2. 其 次 ， 通 过 对 理论 计算 、HTRI 软 件 及 CFD 仿 真 模拟 等 三 种 换 热 器 的 设计 方法 在 熔 盐 换 热 器 的 传 热 性 能 

和 压 降 性 能 方面 进行 的 比较 分 析 ，HTRI 软 件 和 CFD 仿 真 的 温度 计算 结果 与 MSRE 实 验 结果 的 差距 非常 
小 ， 最 大 误差 分 别 为 -1.6% 和 -1.7%。 对 于 壳 侧 压 降 和 管 侧 压 降 ， 理 论 计 算 和 仿真 分 析 的 结果 均 比 实验 结 
果 高 。 

3. 最 后 ， 综 合 本 文 分 析 ，HTRI 软 件 模 拟 和 CFD 念 真 的 结果 与 实验 数据 最 为 接近 ， 而 Kern 方 法 作为 一 种 保 

守 的 理论 计算 方法 ， 与 实验 结果 最 为 偏离 。 其 中 ，HTRI 软 件 相 较 于 CFD 仿 真 更 为 简便 ， 尤 其 是 在 换 热 

# 结 构 是 商业 上 成 熟 的 型 式 时 ，HTRI 软 件 非常 适合 用 作 工 程 设计 参考 。 为 后 续 熔 盐 换 热 器 的 设计 提供 
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